Probing the charge of a quantum dot with a nanomechanical resonator 



H.B. MeerwaldtQG. Labadze, B.H. Schneider, A. Taspinar, Ya.M. Blanter, H.S.J, van der Zant, and G.A. Steel^] 

Kavli Institute of Nanoscience, Delft University of Technology, Lorentzweg 1, 2628 C J Delft, The Netherlands 

(Dated: October 10, 2012) 

We have used the mechanical motion of a carbon nanotube (CNT) as a probe of the average charge 
on a quantum dot. Variations of the resonance frequency and the quality factor are determined by the 
change in average charge on the quantum dot during a mechanical oscillation. The average charge, 
in turn, is influenced by the gate voltage, the bias voltage, and the tunnel rates of the barriers 
to the leads. At bias voltages that exceed the broadening due to tunnel coupling, the resonance 
frequency and quality factor show a double dip as a function of gate voltage. We find that increasing 
the current flowing through the CNT at the Coulomb peak does not increase the damping, but in 
fact decreases damping. Using a model with energy-dependent tunnel rates, we obtain quantitative 
agreement between the experimental observations and the model. We theoretically compare different 
contributions to the single-electron induced nonlinearity, and show that only one term is significant 
for both the Duffing parameter and the mode coupling parameter. We also present additional 
measurements which support the model we develop: Tuning the tunnel barriers of the quantum dot 
to the leads gives a 200-fold decrease of the quality factor. Single-electron tunneling through an 
excited state of the CNT quantum dot also changes the average charge on the quantum dot, bringing 
about a decrease in the resonance frequency. In the Fabry-Perot regime, the absence of charge 
quantization results in a spring behaviour without resonance frequency dips, which could be used, 
for example, to probe the transition from quantized to continuous charge with a nanomechanical 
resonator. 

PACS numbers: 85.85. +j, 81.07.Oj,81.07.De 



I. INTRODUCTION 

Nanomechanical systems^! are studied intensively for 
both their potential applications such as mass sensingp^ 
and for insights into the quantum mechanical ground 
state of a macroscopic objeciP^. Because of their small 
size, nanomechanical resonators are strongly influenced 
by electrostatic forces from single-electron charge effects, 
which allows single-electron transistors to be used as 
sensitive detectors of the deflection of a nanomechani- 
cal beam[LLLL21 demonstrating clear backaction from the 
single-electron forces^. Coupling of these forces to me- 
chanical resonators can also be exploited in mechanical 
single-electron shuttle devices^ to shuttle electrons one- 
by-one through the nanomechanical resonator. 

A carbon nanotube (CNT) is a stiff, bottom-up 
nanomechanical resonator with a large aspect ratio 15 . 
Dissipation in an ultraclean CNT at cryogenic temper- 
atures is low, which results in a high quality factor^ 
and allows investigation into other sources of damping, 
such as nonlineai^ and magnetic damping 18 . At cryo- 
genic temperatures, a quantum dot is formed, embedded 
in the CNT 19 * 20 !, which makes single-electron charge ef- 
fects couple stron gly to the mechanical motion through 
the bending modeP^U. For a CNT quantum dot, the 
effects of damping, spring stiffening and softening, and 
nonlinearity, are completely dominated by single-electron 
charging effects^. 

The interplay between single-electron tunneling and 
mechanical motion has been the topic of many theo- 
retical investigations. The motion of nanomechanical 
resonators is found to have an influence on the elec- 



tron transport t hrough the single-electron transisto j 24 * 25 * 
affecting currentP^ 2 ^ and current noise^HSl Con- 
versely, transport through the single-electron transistor 
by tunneling of single electrons causes backactio n 1 29 * 32 ^ 3 -^ 
on t he mec hanical moti on in the form of frequency 
shifts^ 3 - 5 ^ and damping} 27 1 35 1 37 1 38 1 . Under certain con- 
ditions, single-electron tunneling can lead to a negative 
charging energ}f^ or negative damping, causing instabil- 
ities, where a distinction should be made between the 
low-frequency limit 40 (coo <C T) and the high-frequency 
limit 41 (coq ^> T), relating the mechanical resonance fre- 
quency, cjo, to the single-electron tunnel rate, T. Besides 
causing backaction on the mechanical motion, single- 
electron tunneling is proposed to be used to paramet- 
rically drive the nanomechanical resonator 42 . In the 
Coulomb blockade regime, CNT quantum dot resonators 
in particular are found to have a large electron-vibron 
coupling 43 , and nonlinear restoring forces are found to 
be completely dominated by single-electron tunneling 
effects^El. Further theoretical studies have been per- 
formed on the subject of single-electron shuttlesS3H33 and 
on the coupling between a single-electron transistor and 
a nanomechanical resonator in the quantum regime^H^. 

In this article, we present measurements of single- 
electron effects in CNT resonators in the Coulomb block- 
ade regime and we demonstrate that our experimen- 
tal observations agree quantitatively with the theoreti- 
cal model we develop. This agreement allows the CNT 
nanomechanical resonator to be used as a probe for the 
average charge residing on the CNT quantum dot. Fur- 
thermore, we examine the implications of the established 
model through additional experiments. The layout of the 
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paper is as follows. In section [TTJ the fabrication of the 
ultraclean suspended CNT is described, followed by the 
measurement setup. Section |III| contains the character- 



ization of the CNT device in electrical terms, and the 
influence of single-electron tunneling on the mechanical 
resonance frequency, followed by the novel experimental 
observation of a double frequency dip feature. In sec- 
tion [IV| we develop a model for the dynamics of the me- 
chanical resonator in the presence of Coulomb blockade, 
which explains the presence of such a double frequency 
dip, and make a quantitative comparison to the measured 
data. In section [V| we extend this model to include the 
effect of Coulomb blockade on mechanical damping and, 
again, compare the model to experimental data. In sec- 
tion |VI[ we expand the model further with the descrip- 
tion of single-electron induced nonlinearity, resulting in 
an expression for the Duffing parameter and the mode 
coupling parameter. In section \VU\ we explore the cou- 
pling of the average charge to the mechanical resonator 
by varying the tunnel rates of the quantum dot, by study- 
ing the effects of excited states of the quantum dot, and 
by studying the mechanical resonator in the Fabry-Perot 
regime, in which Coulomb Blockade no longer plays a 
role. All of these additional measurements can be un- 
derstood qualitatively in the context of the model we 
present. 



II. FABRICATION AND MEASUREMENT 
SETUP 



third of the patterned trenches has a CNT across them, 
touching both the source and the drain, thus forming a 
device. Room temperature measurements of current as 
a function of gate voltage are performed for each trench, 
showing semiconductor behaviour for potential devices. 

Figure [TJi shows a schematic diagram of the setup used 
to measure mechanical resonances in the suspended car- 
bon nanotube devices^. The device is mounted at the 
mixing chamber of a 3 He/ 4 He dilution refrigerator with 
a temperature of 20 mK. Filtered twisted-pair cabling is 
used to connect to the source, drain, and gate of the de- 
vice, allowing d.c. voltages to be applied to the source 
and gate. The current flowing through the device is mea- 
sured at the drain. The CNT is driven into motion by 
an a.c. voltage difference between the gate electrode and 
the CNT. This high-frequency signal needed to drive the 
CNT is supplied by an RF source through a coaxial cable. 
At a separation of ~1 cm from the device, the shielding of 
the coaxial cable is removed to form the antenna. We ex- 
pect that the electrostatic coupling between the antenna 
and the nanotube segment itself is much too small to ac- 
tuate the CNT, which was confirmed by a lack of response 
of the CNT quantum dot to a d.c. voltage applied to the 
coax. Instead, the coaxial cable is capacitively coupled 
to the d.c. wires leading to the source, drain, and gate. 
Because of a difference in crosstalk capacitance from the 
coaxial cable to the source, drain, and gate, a.c. volt- 
ages are generated asymmetrically on the three. The a.c. 
voltage difference arising between the gate and the CNT 
then actuates the CNT into motion. 



Ultraclean suspended carbon nanotube devices are fab- 
ricated as follows^. The fabrication begins with a de- 
generately doped silicon wafer with a 285 nm thermal 
oxide. In the first step, the contacts are patterned by 
evaporating 5 nm of tungsten and 25 nm of platinum on 
a patterned double layer of poly(methyl methacrylate) 
(PMMA), and performing lift-off. In the second step, 
a three-layer etch mask, consisting of photoresist, tung- 
sten, and PMMA, is used to etch the trenches between 
the contacts. The trenches are first patterned onto the 
PMMA by electron beam lithography. This pattern is 
transferred onto the tungsten by parallel plate reactive 
ion etching using a mixture of SF6 and helium, where the 
PMMA acts as an etch mask. After this, the photoresist 
is etched by an oxygen plasma, during which the tungsten 
acts as an etch mask. In the third step, trenches between 
the contacts are etched into the silicon oxide, by the same 
mixture of SFq and helium, during which the photoresist 
acts as an etch mask. To improve wire-bonding, a layer 
of 10 nm of chromium and 80 nm of platinum is evapo- 
rated onto the bondpads, followed by a sputtered layer 
of 20 nm of silicon. In the final step, catalyst islands 52 
are deposited onto holes patterned in a double layer of 
PMMA. The sample is now placed in a CVD oven, where 
CNTs grow out of the catalyst particles in a mixture of 
hydrogen and methane at a temperature of 900 °C. As 
the CNTs grow in a random direction, approximately one 



III. ELECTRICAL AND MECHANICAL 
CHARACTERISTICS 

Figure [TJd shows the current through a small-bandgap 
CNT with a suspended length of 600 nm as a function of 
gate voltage. For gate voltages below V g = 0.4 V, the de- 
vice is doped with holes and weak scattering at the metal- 
CNT interface at the edge of the trench results in conduc- 
tance that is modulated by Fabry-Perot interferences^ 
of the hole wavefunction. For 0.4 < V g < 0.5, the 
Fermi level lies in the bandgap and the current is sup- 
pressed. From the distance in gate voltage from electron 
and hole conduction, together with the coupling factor 
a = Cg/Ctot = 0.38 determined from the Coulomb dia- 
monds, we estimate the bandgap to be E gap = 58 meV. 
Above V g = 0.5 V, Coulomb oscillations are visible. Now, 
electrons tunnel onto the CNT through tunnel barriers, 
which arise from the p-n junctions between the segments 
of the CNT near the W/Pt metal near the edge of the 
trench and the CNT. The small capacitance between the 
quantum dot and the three terminals gives rise to a large 
charging energy of Ec = 9.6 meV. For V g > 0.7 V, the 
Coulomb peaks are increasingly more smeared out. For 
larger electron doping, the p-n junctions become nar- 
rower and the tunnel barriers to the quantum dot become 
more transparent 51 . 



We measure the mechanical bending mode resonances 
of the CNT by actuating it into motion and measuring 
the d.c. current. When the drive frequency matches 
the resonance frequency, the CNT resonates, causing the 
capacitance between the CNT and the gate to oscillate 
with a large amplitude. This oscillating capacitance ef- 
fectively induces an oscillating gate voltage. The nonlin- 
earity of the Coulomb peak allows the oscillating effective 
gate voltage from the motion to be rectified into a d.c. 
current!^. 

As reported earlier, the mechanical resonance fre- 
quency of the CNT is strongly influenced by single- 
electron tunneling^. The middle panel of Fig. [lp shows 
the change in current due to the mechanical motion, A/, 
as a function of the drive frequency and gate voltage on 
a Coulomb peak at Vb = 0.2 mV. On a Coulomb peak, 
single-electron tunneling leads to a dip in the resonance 
frequency. In the top panel of Fig. a linecut is shown 
of current versus gate voltage at a drive frequency chosen 
to be far from the resonance frequency. In the right panel 
of Fig. [l£, a linecut is shown of the current versus the 
drive frequency at one gate voltage, denoted by the black 
dashed line in the middle panel, showing the mechanical 
resonance as a change in the d.c. current. 

In order to qualitatively understand these dips in fre- 
quency, three essential elements are required. The first 
is that the motion of the CNT changes the charge on 
the quantum dot, which in turn changes the electro- 
static force on the CNT. Thus, the CNT experiences a 
displacement-dependent electrostatic force. Because the 
electron tunnel rates T are much faster than the mechan- 
ical frequency f (T ~ 450 GHz, f ~ 300 MHz), the 
mechanical motion sees a displacement-dependent force 
from the charge averaged over many tunnel events. This 
displacement-dependent force induces a reduction of the 
spring constant, which changes the mechanical frequency. 

The sign of the frequency shifts can be understood by 
realizing that electrostatic forces do not act as a restoring 
force, but instead as an anti-restoring force: if the CNT 
is pulled away from its equilibrium position towards the 
gate, for example, the electrostatic force will increase and 
tend to pull the CNT with more force towards the gate. 
Thus, electrostatic forces generally result in a decrease of 
the net spring constant. 

Finally, to understand the gate voltage dependence, it 
is illustrative to examine the differential capacitance of 
the quantum dot to the gate. Due to the quantized charge 
on the island of the dot, the quantum dot shows zero dif- 
ferential capacitance when it is in the Coulomb valleys: 
the charge is fixed independent of gate voltage. At the 
position of the Coulomb peaks, the charge on the island 
undergoes a discrete step of one electron. In the absence 
of tunnel or temperature broadening of this transition, 
this step would be infinitely sharp, and the differential 
capacitance would diverge. The electrostatic spring ef- 
fects discussed in the previous paragraph are determined 
by the differential capacitances for small amplitudes of 
motion, and the diverging differential capacitance of the 
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FIG. 1. (a) Schematic drawing of the chip geometry, antenna, 
and measurement electronics. The CNT acts as a doubly- 
clamped beam resonator with a displacement x and is driven 
because of an asymmetrical capacitive coupling of the radio 
frequency coaxial cable to the source, drain, and gate, (b) 
Current versus gate voltage at Vb = 0.3 mV showing, for in- 
creasing gate voltage, Fabry-Perot oscillations, then a small 
bandgap, and finally Coulomb oscillations with the increasing 
tunnel coupling to the leads opening up the quantum dot as 
the gate voltage is increased further, (c) Top panel: Current 
as a function of gate voltage, showing a Coulomb peak. Mid- 
dle panel: Change in current as a function of drive frequency 
and gate voltage, showing a dip in resonance frequency across 
the Coulomb peak. Right panel: Current as a function of 
drive frequency at a gate voltage, V g = 508.5 mV, as denoted 
by the black dashed line in the middle panel, showing the me- 
chanical resonance as a decrease in current around / = 287.33 
MHz. 



quantum dot at the positions of the Coulomb peaks leads 
to a diverging softening renormalization of the net me- 
chanical spring constant, resulting in the dips in mechan- 
ical frequency observed in the data. 

In this article, we explore in detail how single-electron 
tunneling affects the mechanical motion of the CNT and 
develop a quantitative model to describe our results. In 
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contrast to earlier work^, here we map out the behaviour 
of the mechanical resonator for bias and gate voltages 
covering the single-electron tunneling region between two 
charge states. The stability diagram in Fig. [2^i shows the 
differential conductance, dl/dVb, as a function of gate 
voltage and bias voltage for the transition from 1 to 2 
electrons. Coulomb blockade is visible in blue, whereas 
single-electron tunneling takes place in the red and white 
regions. 

At a low bias voltage, Fig. [2}d shows, in blue, the ex- 
perimentally obtained normalized current, \(I — Iq)\/\I — 
Io\max, as a function of gate voltage and drive frequency, 
where Io is the current off mechanical resonance. The 
bias voltage of 0.17 mV for this figure is denoted by the 
upper white dashed line in Fig. [2^i. The mechanical reso- 
nance is visible as an increase in the normalized current. 
Around V g offset = 1 mV, where V g offset = V g - 0.565 
V, the mechanical resonance frequency shows a dip, as 
was demonstrated in previous measurements. 

When the bias voltage is increased, the resonance fre- 
quency of the CNT exhibits substantially different be- 
haviour. Figure^ shows the measured normalized cur- 
rent as a function of gate voltage and drive frequency as a 
color plot at a bias voltage of —0.55 mV, corresponding to 
the lower white dashed line in Fig. Instead of one dip, 
two dips in resonance frequency are visible. The presence 
of such a peculiar double dip feature in the mechanical 
resonance frequency forms the motivation for the work 
in the following sections, in which we establish a quanti- 
tative model for the mechanical resonance frequency and 
quality factor in the presence of Coulomb blockade, and 
further explore the coupling of the mechanical resonator 
dynamics to the average charge of the quantum dot. 



IV. MODELING SINGLE-ELECTRON SPRING 
EFFECTS 

This section describes a model developed for the cou- 
pling of the mechanical resonator to the Coulomb block- 
aded quantum dot, and the quantitative fitting of the 
experimentally observed frequency dips. In subsection 
|IV A[ a model for the average charge on the CNT is de- 
scribed. In subsection |IVB[ we derive a model for how 
the average charge leads to a displacement-dependent 
force that causes the softening of the CNT spring. In 
subsection |IVC[ fits of the experimentally obtained res- 
onance frequency and current are performed using the 
established model. 



A. Model for the average charge on a quantum dot 

The average charge on the CNT at a charge transition 
can be determined as follows. Single-electron tunneling 
onto or off the CNT is a stochastic process, where the 
amount of time an electron spends on the source or the 
CNT is determined, respectively, by the tunneling rates 
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FIG. 2. (a) Stability diagram: dl/dVb as a function of bias 
and gate voltage for the charge transition from 1 to 2 elec- 
trons. The linecuts, at which Figs. [2Jd and c are taken, 
are indicated by the white dashed lines. The black dashed 
line denotes the gate voltage at which the middle panel 
of Fig. Et is taken, (b, c) Measured normalized current, 
|(7 — Jo)17l^ — Io\max, as a function of drive frequency and 
gate voltage offset at (b) V& = 0.17 mV, showing a decrease 
in resonance frequency around Vg )0 ffset = 1 mV, and at (c) 
Vb = —0.55 mV, showing a double dip in resonance frequency. 

T + and r~. The average occupation, (N), of charges on 
the CNT for the charge transition from Nq to Nq + 1 
electrons is given by: 

(N)=N + ^, (1) 

1 tot 

where T tot = T+ + T~ . 
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The tunneling rates, determining the average occupa- 
tion, are modeled as follows. Traditionally, tunneling 
onto and off a quantum dot is describe d 54 * 55 * using energy- 
independent tunnel rates. The average occupation and 
the current are determined by using energy-independent 
tunnel rates and the overlap of the density of available 
states of the quantum dot and the leads. Such an ap- 
proximation of energy-independent tunnel rates is valid 
when the tunnel barriers are sufficiently high, so that a 
change in either bias or gate voltage does not cause a 
significant change in the barrier height or width. 

In this work instead, we use energy-dependent tunnel 
rates, assuming that the tunnel barriers are not high. 
This assumption is supported by the strong dependence 
of the tunnel barriers on the gate voltage (cf Fig. [I}}). 
The influence of bias and gate voltage on the tunnel rates 
is demonstrated in Figs. [3^l and b. At zero bias voltage, 
Fig. [3^i shows the p-n junctions as tunnel barriers be- 
tween the CNT and the leads. When the bias voltage 
is increased, Fig. [3}d illustrates how the higher electro- 
static potential of the source increases the tunnel barrier 
at the source, leading to a lower tunnel rate, whereas 
the tunnel barrier at the drain is decreased, resulting in 
a higher tunnel rate. In our device, the tunnel barriers 
are not formed by steps in the potential but instead by 
p-n junctions in the CNT. Figure ^ shows how tunnel 
barriers arise from p-n junctions formed in the CNT at 
the metal-CNT interface, as the chemical potential of the 
CNT is decreased by a positive gate voltage. The height 
of the tunnel barriers is determined by the bandgap of 
the CNT, E gap . 

To calculate the tunnel rates, the density of available 
states of the leads and the quantum dot are used, taking 
the energy dependence of the tunnel rates into account. 
The density of available states of the level in the CNT 
quantum dot, D(/i), caused by tunnel coupling to the 
leads, is described with a Lorentzian lineshape. 

= 2^(^-Mcivr) 2 + ^r 6 w/2) 2 - (2) 

Here, lie NT is the chemical potential of the CNT. The 
broadening T^road gives the full width at half maximum. 
At T = 20 mK, we are in the regime where TiViy roa d ^> 
ksT. The density of available states of the left and right 
lead, /l,r(E)^ is now modeled by a step function: 

>»<*>-{ J; i^SE, (3) 

where E is the energy of the electron, and Hl,r is the 
chemical potential of the left (right) lead. The tunnel 
rates are determined by the overlap of the density of 
available states of the CNT and the leads, and by their 
relative chemical potential: 

r*. = G + ; — (^))- 

(4) 



(b) 




(c) 
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FIG. 3. (a) Energy diagram at zero bias voltage, with the p-n 
junctions represented by square tunnel barriers between the 
CNT and the leads, (b) Energy diagram at finite bias voltage, 
showing an increase of the height of the tunnel barrier to the 
source relative to the energy of the tunneling electron, and a 
decrease in the height of the barrier to the drain, (c) Diagram 
of the CNT band structure with two p-n junctions induced 
through a positive gate voltage. 



Here, the energy-dependent tunneling is reflected by the 
exponential factor, which is shown theoretically 56 to oc- 
cur for tunnel barriers brought about by an electrostatic 
potential, where the barrier height is low and the bar- 
rier width is large. The parameters cll,r depend on the 
height and width of the tunnel barriers. The parameters 
bL,R describe the triangular profile that the tunnel bar- 
riers have, which is there at zero bias voltage (cf Fig. [3^) 
and is additionally changed by altering the bias and gate 
voltage (cf Fig. [3)3). The difference in chemical potential 
between the CNT and the left or right lead is denoted by 
AfiL,R = He nt — Hl,r- The last factor in Eq. [4] arises 
from the broadening due to tunnel coupling and is deter- 
mined from the overlap of the density of available states 
of the left or right electrode with that of the CNT. 

The calculated average occupation, at a bias voltage 
of Vb = 0.17 mV, corresponding to the white dashed line 
labeled (b) in Fig. [2^, is shown in the top panel of Fig. 
|4|l as a function of gate voltage at the charge transition 
from 1 to 2 electrons. At zero broadening, T^road = 0, 
the green line shows the average occupation increase in 
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two discrete steps. The left inset of Fig. [4^i shows an en- 
ergy diagram illustrating zero broadening. The three dis- 
crete plateaus in average charge correspond to Coulomb 
blockade, single-electron tunneling, and Coulomb block- 
ade again. At a broadening of T^road = 450 GHz, for 
which hT b roa( i > e\Vb\, the red line shows the average 
occupation increases monotonically in a single step. The 
right inset of Fig. [4^i illustrates how the double step 
in average occupation, seen at zero broadening, is com- 
pletely smeared out. Here, e is the elementary charge, ks 
is Boltzmann's constant, and Ti is the reduced Planck's 
constant. 

At a bias voltage of V5 = —0.55 mV, corresponding to 
the white dashed line labeled (c) in Fig. [2^1, the calcu- 
lated average occupation is shown in the top panel of Fig. 
[4}d as a function of gate voltage. Both at a broadening of 
Abroad = (green) and at a broadening of Y^road — 450 
GHz (red), two steps and three plateaus are visible. The 
insets of the top panel of Fig. [4}} show how all allowed 
states on the CNT are inside the bias window, both for 
zero broadening (left) and finite broadening (right). Con- 
trary to low bias, the two steps are not smeared into a 
single step at finite broadening. 

For different bias voltages, the calculated average occu- 
pation is shown in the top panel of Fig. as a function 
of gate voltage for different bias voltages. At low bias, 
the charge transition takes place in a small range in gate 
voltage. This gives rise to a large slope for the aver- 
age occupation with respect to gate voltage. As the bias 
voltage is increased, the range in gate voltage increases. 
When the bias voltage is larger than the broadening, the 
average occupation increases in two steps. 



B. Model for the single-electron spring 

Using the average charge as determined in the previous 
subsection, the shift in resonance frequency is modeled 
as follows. The electrostatic force acting on the CNT 
depends on the voltage difference between the CNT and 
the gate electrode 57 : 



F c 



NT 



l dCg 

2~dx 



(V 9 



Vcnt) 2 \ 



(5) 



where C g is the capacitance between the CNT and the 
gate, x is the displacement of the fundamental mode of 
the CNT, and V g and Vcnt are the voltages on the gate 
and the CNT, respectively. The voltage on the CNT is 
determined by the control charge, q c = C g V g + CsVs + 
CdVd, which is the charge that would be on the CNT 
in the absence of Coulomb blockade, and the average 
occupation, (TV), of charges residing on the CNT: 



Vcnt 



Qc ~ e(N) 

C+ot 



(6) 



Here, Cs,d an d Vs,d are the capacitances to and the 
voltages on the source and drain respectively, and Ctot — 



C g + Cs + Cd- I n our case, the bias voltage is applied to 
the source, Vs = V5, and the drain electrode is grounded, 
Vd — 0. Because the charge N on the dot fluctuates 
stochastically between No and 7Vo + l at a rate T to t ^> /o? 
the mechanical motion experiences a voltage on the CNT 
due to an average occupation, (N). 

The spring constant and, consequently, the resonance 
frequency of the CNT are determined by the change in 
force acting on the CNT per unit displacement: 



2racJoAcJo = Afc : 



dF 

dx ' 



(7) 



where m is the mass of the CNT, and Aujq <C ojo- The full 
derivative of the force with respect to the displacement of 
the CNT is expanded into partial derivatives with respect 
to displacement and gate voltage: 



dF(x,V g 
dx 



dF(x,V g ) _ dC g dq c dV g dF(x,V g 
dx 

ld 2 C a 



(V 9 



2 dx 2 
V 9 (V 9 -Vcnt) 



dx dC g dq c 
Vcnt) 2 



dV a 



dCgV d(Vg-V CN T) 



dx 



dV n 



(8) 



The first term gives rise to a softening spring effect due 
to the capacitive force between the CNT and the gate. 
The second term takes into account the influence of the 
displacement on the control charge through the gate ca- 
pacitance. In turn, the influence of the control charge 
on the force is incorporated through the gate voltage. In 
this article, we focus on changes in the spring constant 
that occur rapidly with gate voltage and we show that 
the experimental features can be captured using only the 
second term. 

Combining Eq. (|6| and the second term of leads 
to the following expression for the change in resonance 
frequency due to a changing average charged 



2muJoCt ( 



Vcnt) ( dC g 
dx 



Ctot 

~c7 



e d(N) 
~C a dV a 



( 9 ) 

Because of the rightmost minus sign and the fact that 
d (TV) jdVg > 0, the changing average charge on the CNT 
leads to a softening spring effect. The resonance fre- 
quency of the CNT decreases more when the mechanical 
oscillation causes a larger change in average charge, ex- 
pressed by d (N) jdVg. 



C. Fitting of the experimental resonance frequency 
shift and discussion 



To verify the model established in the previous two 
paragraphs, we perform a quantitative fit on the experi- 
mental data shown earlier in Fig. [2] The fits are accom- 
plished by using Eq. [9] and a single set of parameters for 
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FIG. 4. (a, top) Calculated average occupation, (N), as a function of gate voltage at V& = 0.17 mV, showing two discrete 
steps for Tbroad — (green line), and a smeared-out monotonic increase for Tbroad — 450 GHz (red line). Insets: energy 
diagrams at low bias for zero broadening (left) and finite broadening (right), (a, b, c, middle) Measured normalized current, 
\{I — Io) | /| J — Jo | max, as a function of drive frequency and (a, blue) gate voltage offset at V& = 0.17 mV (b, blue) gate voltage 
offset at Vb — —0.55 mV, and (c, red) bias voltage at a gate voltage offset V g , ffset = 0.04 mV and the corresponding calculated 
resonance frequency (red line, red line, blue line), showing (a) a decrease in resonance frequency around V g offset = 0, (b) a 
double dip in resonance frequency, (c) a decrease in resonance frequency for low bias voltages, (a, b, c, bottom) Measured 
(blue dots) and calculated (red line) current as a function of gate voltage at (a) Vb = 0.17 mV, (b) Vb = —0.55 mV, and (c) as 
a function of bias voltage at V g , ffset = 0.04 mV. (b, top) Calculated average occupation, (N), as a function of gate voltage 
at Vb = —0.55 mV, showing two discrete steps for both Tbroad = (green line) and for Tbroad = 450 GHz (red line). Insets: 
energy diagrams at high bias for zero broadening (left) and finite broadening (right), (c, top) Calculated average occupation, 
(iV), as a function of gate voltage and bias voltage at Tbroad = 450 GHz, showing the average occupation, (iV), from increasing 
in a single step at zero bias to increasing in two steps at Vb — 1 mV. 



all figures. The values chosen for the parameters can be 
found in appendix [A! A value for d (N) / dV g is obtained 
by numerically differentiating the calculated average oc- 
cupation, as is displayed in the top panel of Fig. [4^i, with 
respect to gate voltage. 

The red graphs in the middle panel of Fig. [4ji and 
b show the calculated change in resonance frequency as 
a function of gate voltage. At Vb = 0.17 mV, Fig. [4^i 
shows quantitative agreement between the measurement 
and the model describing the frequency dips. The fre- 
quency dip corresponds to the largest slope of the aver- 
age charge with gate voltage. At Vb — —0.55 mV, Fig. 
[4}d demonstrates how the model reproduces the experi- 
mentally observed double dip structure. At the plateau 
in average charge, the mechanical oscillation only brings 
about a small change in average charge, and the reso- 
nance frequency returns towards its original value. 



To investigate the reduction in the resonance frequency 
as a function of bias voltage, a vertical linecut is taken, 
at a constant gate voltage through the charge degeneracy 
point in the stability diagram. The middle panel of Fig. 

shows the measured current, but now as a function 
of bias voltage, not gate voltage, and drive frequency. 
The calculated resonance frequency in blue shows excel- 
lent quantitative agreement with the measurement. At 
low bias voltages, the narrow charge transition leads to 
a large change in average charge due to the mechani- 
cal oscillation and a large change in the displacement- 
dependent force, resulting in a large decrease in resonance 
frequency. At higher, both positive and negative, bias 
voltage, the slope in average occupation, d(N)/dV g , be- 
comes less, as the plateau arises. This leads to a smaller 
change in the displacement-dependent force and a smaller 
decrease in the resonance frequency. 
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The experimentally obtained current flowing through 
the CNT is shown with blue dots in the bottom panels 
of Fig. gat V b = 0.17 mV (a), V b = -0.55 mV (b), and 
as a function of bias voltage (c). For Vb = 0.17 mV, a 
Coulomb peak is visible. For Vb = —0.55 mV, the plateau 
in current shows that the bias voltage, V b , is larger than 
the broadening, T broad- As a function of bias, the bottom 
panel of Fig. shows no Coulomb blockade, since the 
vertical linecut in bias voltage exactly passes through the 
charge degeneracy point. The absence of a saturation of 
the current at high bias voltage indicates low tunnel bar- 
rier heights and supports the choice of energy-dependent 
tunnel rates. As an independent examination of the cho- 
sen tunnel rates, the measured current is fitted using the 
same parameters as for the frequency dips (red lines). 
The current is calculated using: 



r£rv 



(10) 



For intermediate bias voltages, Fig. [5] illustrates the 
transition from the double dip structure to a single dip. 
Using the same set of parameters as for Fig. [4j the calcu- 
lated resonance frequency follows the measurement well. 
The quantitative agreement between the experiment and 
the calculation of the resonance frequency and the cur- 
rent for different bias voltages demonstrates the consis- 
tency of the single set of parameters used for the calcu- 
lation. 



V. SINGLE-ELECTRON DAMPING 

In this section, the influence of the changing average 
charge on the mechanical quality factor of the CNT is 
investigated. As with the resonance frequency in the 
previous section, the quality factor is determined experi- 
mentally and by using a model, for low bias voltage, high 
bias voltage, and as a function of bias voltage. To show 
the consistency of the used parameters, the correspond- 
ing resonance frequency is also displayed. 

Figure [6ji shows the stability diagram of the charge 
transition from 7 to 8 electrons at a parallel magnetic 
field of 8 T. The current flowing through the CNT is 
measured as a function of bias and gate voltage, and 
the derivative dl/dVb is determined numerically. The 
letters b and d denote, respectively, the low and high 
bias voltage, at which Figs. [BJd and|6]i were taken. The 
letter f denotes the gate voltage at which Fig. [6]f was 
taken. 



A. Experimental observation of single-electron 
damping 

The experimentally determined mechanical quality fac- 
tor at a low bias voltage of V b = 0.3 mV is displayed in 
blue dots in the top panel of Fig. [6]o. The quality factor is 
obtained by measuring the current versus drive frequency 



10 times and averaging. The drive power is adjusted for 
each quality factor measurement such that the frequency 
response shows a lineshape which is as close as possible 
to Lorentzian, but which is not obscured by noise. By 
working at low powers in the linear response regime, we 
minimize the likelihood that non-linear damping termd^ 
play a significant role. The change in current, A/, due to 
the mechanical motion in our detection scheme is propor- 
tional to the amplitude of the displacement squared. At 
the sides of the Coulomb peak, the lack of nonlinearity of 
the current with gate voltage prevents the measurement 
of the mechanical motion through a change in current, 
and consequently a determination of the quality factor. 
A fit with the following Lorentzian function is performed 
to extract the quality factor. 



AJ(cj) = A- 



(("o 



,2\2 



(11) 



Here, Q is the mechanical quality factor, uj is the drive 

2 



frequency, and A = \ J W^ x<2 * s a ^ parameter, 

which incorporates the amplitude of oscillation x and the 
electromechanical coupling. At a low bias of Vb — 0.3 
mV, a large decrease in the quality factor of two orders 
of magnitude is visible, spanning several millivolts of gate 
voltage corresponding to single-electron tunneling. The 
top panel of Fig. [6]i shows the experimentally obtained 
quality factor in blue dots at a high bias voltage of Vb = 
1.1 mV. Here, similarly to the double frequency dip at 
high bias in the middle panel of Fig. [4|d, two dips in 
quality factor are visible. 

To investigate the relation between current and single- 
electron damping of the mechanical motion, Fig. [6^ 
shows the current as a function of bias voltage. The 
gate voltage is chosen such that the graph is a vertical 
linecut through the charge degeneracy point. The top 
panel of Fig. [6]f shows the corresponding quality factor 
as a function of bias voltage, again at a gate voltage of 
V g = 0.9002 V. The quality factor exhibits a dip at low 
bias voltages and returns to its original value at higher, 
both negative and positive, bias voltage. We emphasize 
that the highest current flowing through the CNT does 
not correspond to the largest reduction in quality factor. 
Instead, in fact, when the gate voltage is positioned at a 
Coulomb peak, the quality factor actually goes up when 
increasing the current through the device. 

Figure [B^ shows the basic concept of the damping 
mechanism, illustrated with zero broadening and one 
electron tunneling event during a mechanical oscillation. 
Let us first consider the zero bias voltage case. The me- 
chanical motion brings the chemical potential of the CNT 
below and above the Fermi level of the leads. The asym- 
metry of the Fermi sea in the leads implies that electrons 
can only tunnel onto the CNT when the level of the CNT 
is below the Fermi energy of the leads, and can only tun- 
nel off when it is above. As the mechanical oscillation 
pushes the level downwards through the Fermi energy, 
there is a small retardation in the time when the elec- 
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FIG. 5. Measured normalized current (blue) and calculated resonance frequency (red line) as a function of gate voltage offset 
for different bias voltages, showing the transition of the double dip structure into a single dip as the bias voltage becomes less 
negative. 



tron tunnels on, given by 1/T. As the level is pushed 
back upwards in electrostatic energy by the motion, the 
electron tunnels off, again with a small time delay. The 
net result is that electrons are pumped from below the 
Fermi energy to above it, extracting energy from the me- 
chanical motion, and thus resulting in mechanical damp- 
ing. Although the retardation time is small compared to 
the mechanical frequency (r ^> cjq), the damping is still 
large due to the large electro-mechanical coupling. At 
finite bias voltage, the same picture applies, but now the 
damping occurs when the level of the CNT passes by the 
Fermi energy of each of the left and right lead separately. 

At zero temperature and with no broadening of the 
transition from the tunnel coupling to the leads, the dip 
in quality factor would be infinitely narrow. Including 
finite temperature and tunnel coupling of the level to the 
leads, the dip in quality factor acquires a finite width, 
and is proportional to the change in average occupation 
with gate voltage, d(N) / 'dV g , as shown in Eq. 
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the next subsection. The fact that the quality factor is 
determined in this way also explains the observed sin- 
gle and double dip structures in the quality factor. At 
low bias voltage, the charge transition takes place in a 
single step and consequently there is a single dip in the 



quality factor. As the bias voltage is increased, the single 
step changes into two steps with a plateau between them. 
On the plateau, the change in average occupation with 
gate voltage, d(N)/dV g , is smaller, leading to a smaller 
retarded single-electron force acting on the CNT and a 
smaller reduction of the quality factor. Although a larger 
current is due to more tunnel events, the CNT does not 
perform extra work at high bias. As more electrons pass 
the CNT and reside on it, the retardation time is re- 
duced. What causes the reduction in quality factor is the 
pumping of the electrons from below the Fermi energy to 
above it, and the retardation between the electrostatic 
force and the mechanical motion. 



B. Model for single-electron damping 

Using the following model, the quality factor is calcu- 
lated, as shown in the red lines in the top panels of Figs. 
[6]o, [6]l, and(6f. We look at the limit where there are many 
tunneling events per mechanical oscillation, Y tot ^> / . 
As derived in Refs. [27] and SHI the total damping has 
an intrinsic contribution and a contribution due to the 
displacement-dependent force associated with tunneling 
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electrons: 



UJ 



F stoch^q 



1 dC g d(N} 



mC a 



■ tot 



dx dV a 



(12) 



where Q to t is the total quality factor and Qi n t is the 
intrinsic quality factor in the absence of tunneling elec- 
trons. The stochastic force experienced by the CNT, 
Fstoch = F(No + l) — F(No), is the difference between the 
force experienced at Nq and Nq + 1 electrons. Assuming 
that both source and drain voltage are much smaller than 
the d.c. gate voltage, V 8 ,Vd <C V g dc : , and taking only into 
account the d.c. force acting on the electrons as the CNT 
oscillates (which is valid as long as Vg C <C V g dc ), Ref. [271 
gives the following expression for the stochastic force: 



1 dC a 



stock 



^tot 



dx 



(2e(C s + C D )V g + e 2 (2N + l)). 



(13) 

Equation (12) illustrates how single-electron damping is 
increased by a large change in average occupation be- 
cause of mechanical motion, d(N)/dV g , but is reduced 
by a large total tunnel rate, T totl as the retardation time 
decreases. In the case that exponential tunnel rates cause 
d(N)/dV g to become negative, single-electron tunneling 
pumps the mechanical motion, leading to self-sustained 
oscillation^. 



Fitting of the experimental quality factor and 
discussion 



The calculated and the measured quality factors, at 
low bias, are compared in the top panel of Fig. [6)3, 
showing quantitative agreement below V g offset = 1 mV, 
where V gj0 ff set = V g — 0.9 V. At high bias, the top panel 
of Fig. |6]i shows an excellent quantitative agreement be- 
tween the model and the measurement. As a function 
of bias in the top panel of Fig. [6]f, the calculated qual- 
ity factors match the measured quality factors at bias 
voltages below 1 mV. Above 1 mV, an excited state, as 
evident in the current in Fig. [6^, may be the cause of 
extra damping. 

The bottom panels of Figs. [6}d, [6]i, and [6]f show, in 
blue dots, the resonance frequency extracted from the 
frequency responses alongside the quality factors. The 
decrease in quality factor is accompanied by a decrease 
in the resonance frequency. The red lines show the reso- 
nance frequency calculated with the model as described 
in the previous section, using the same tunnel rates as 
for the quality factors. Without changing any fit pa- 
rameters, we have qualitative agreement as the quality 
factor dips and the resonance frequency dips occur at 
corresponding gate voltages. It is not fully understood 
why agreement between the calculated and measured res- 
onance frequency is not quantitative. 



VI. NONLINEAR RESTORING FORCES DUE 
TO SINGLE-ELECTRON TUNNELING 

In the previous sections, we described linear correc- 
tions to the restoring force due to single-electron tunnel- 
ing, which change the spring constant and cause damp- 
ing. In this section, we cover nonlinear corrections due to 
single-electron charge effects^, resulting in expressions 
for the Duffing parameter and the mode coupling pa- 
rameter. 



A. Duffing nonlinearity due to single-electron 
tunneling 

Adding nonlinear corrections, the single-electron force 
is expanded with respect to the displacement, x, as fol- 
lows: 

Fcnt = -AksETX - Psetx 2 - asETX 3 - (14) 

The Pset parameter and the Duffing parameter, qlset-, 
are calculated as derivatives of the single-electron force: 



PsET 



asET 



1 d 2 FcNT 

"2 dx 2 

1 d 3 FcNT 

6 dx 3 



(15) 
(16) 



Expanding the full derivatives into partial partial deriva- 
tives as d/dx — d / dx+(V g /C g )(dC g / dx)d / dV g , we arrive 
at: 



_ 1 V g dC g 1 (V g dC g 

Pset - - - 2 F XX - cT^F xg - - I -— ) F gg 



(17) 



„ -- l -F - l -Y^ d Sk F 

L^SET — n r xxx ri 7 r xxq 

6 2 C g dx 

2 



2 {(■„ dx J 6 \C g dx ) 

(18) 

Here, the subscripts of F denote differentiation with re- 
spect to displacement, x, or gate voltage, g. 

The quadratic nonlinearity, ft set, renormalizes^ 
asET, leading to an effective Duffing parameter, 

a e ff,SET- 



10 (3 2 

&eff,SET = ®SET + &OLR.SET = &SET 7T ~ o 



SET 



■ (19) 



For the parameters of the device we study here, there is 
a leading order dominant contribution, oP SETl given by: 



a 



1 fdC g \ V V g Y e(V g - Vcnt) d 3 (N) 



SET 



6 V dx J \C, 



(20) 
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FIG. 6. (a) Stability diagram at a parallel magnetic field of 8 T and V g — 0.9 V, showing the charge transition from 7 to 
8 electrons, with the gate and bias voltages at which Figs. [SJd, |SJl, and [§J were taken, indicated by black and blue dashed 
lines, (b, d, f) Comparison between measurement (blue dots) and model (red lines) of the quality factor (top) and the shift in 
resonance frequency (bottom) as a function of gate voltage (b, d) and bias voltage (f), showing, (b) at V& = 0.3 a single dip, 
(d) at Vb = 1-1 mV, a double dip, and (f) at V g offset = mV, a large decrease in quality factor and resonance frequency at 
low bias voltages, (c) Energy diagram (exaggerated) illustrating the damping mechanism, showing the asymmetry of electrons 
tunneling onto the CNT at low chemical potential and off the CNT at high chemical potential, (e) Measured (blue dots) and 
calculated (red line) current as a function of bias voltage at V g , Q ffset = V. 



which is one of the terms arising from 

-l/6((Vg/Cg)(dCg/dx)) S Fggg ill Eq. M U Slllg 

the same set of parameters as for Figs. ]4]~and [5j we 
repeat the calculated average occupation as displayed 
in Fig. |4^i, and we plot ot e ff,SET and oP SET together in 
Fig.JTl The figure shows that the other contributions in 
Eq. [l8|and the renormalization, Aap^sET, due to ft set-, 
play no significant role in our device. For completeness, 
the other contributions, which are more than one order 
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of magnitude smaller than ql% et , are plotted in Fig 
in appendix [B] The switching of the sign of the Duffing 
parameter a e ff y sET, clue to the third derivative of the 
average occupation with respect to gate voltage, is 
visible in Fig. [7] and has been observed experimentally 
previously^. 



B. Mode coupling due to single-electron tunneling 

Recentl}^, different bending modes in CNT quantum 
dot resonators have experimentally been shown to have 
strong coupling, resulting in a shift in the resonance fre- 
quency of one mode due to the resonance of another 
mode. In contrast to top-down micromechanical beams^Sl 
and CNTs outside the Coulomb blockade regime 61 , mode 
coupling in CNT quantum dot resonators is not domi- 
nated by tension, but, instead, by single-electron charge 



effects. In this section, we establish a theoretical frame- 
work for single-electron mode coupling. We focus on 
the leading order nonlinear contribution to the restor- 



ing force, arising from Fg 



SET 



i SET . 



c 3 , and write it in 
terms of the change in capacitance, 5C 9 , due to mechan- 
ical motion: 



r SET 



-B(V g 



d 3 (N)dC g 
dV g s dx 



(21) 



where, for brevity, we capture the slowly varying depen- 
dence on gate voltage in B(V g ) = \ (^) ^clT^ ' 
The change in capacitance due to mechanical motion de- 
pends on the mode shape: 



SC g = J ^ l^u n (z)cos(u n t)jdz, 



(22) 



where dc g /du is the change in capacitance per unit length 
with respect to the amplitude u{z) at position z along 
the CNT. We use the infinite cylinder parallel to a plate 
model 6 -, and assume that dc g /du is independent of the 
position, z, along the CNT. The position-dependent am- 
plitude, u n {z) is given by: 



(23) 



where the mode shapes, £ n (z), are orthonormalized as 
Jo £i( z )£j( z )dz = LSij with Sij the Kronecker delta, such 
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FIG. 7. (top) Average occupation, (iV), as a function of gate 
voltage, identical to the top panel of Fig. (bottom) Cal- 
culated nonlinearity parameter as a function of gate voltage 
at Vb = 0.17 mV across the Coulomb peak depicted in Fig. 
[2] with its sign nipping from negative to positive to nega- 
tive again, showing the significance of the a° SET contribution 
in a e ff,sET- Other contributions, including the renormaliza- 
tion, Aotp } sET, due to /3set, are shown to be negligible. 



that the amplitude x represents the root-mean-squared 
amplitude over the length of the CNT (not over time). 
With this definition of the amplitude^, the mass of each 
mode is given by the total mass and the spring constant 
is given by k n = muj^. The change in gate capacitance 
due to the mechanical motion is simplified to: 



SC g = L^- y^x n q n cos(u n t), 



(24) 



with the first four elements of a n = (1/L) J ^ n (z)dz 
calculated^ for a bending-rigidity dominated beam to be 
0.83, 0, 0.36, and 0. The leading order nonlinear term, 
Fg ET is expanded as: 

F$ ET = -B(V g )?^C({x i }) 



dV* 



2 2 

T 



+ ]-cos(2u! m t) I a n x n cos w n t 



(25) 



where 



have 



used 



C({xi}) 



(Ldc g /du) A (Y^i a i X i)/ Z) x i' 

We address several terms in this expression. The terms 
containing (1/2) cos(2cj m £) give rise to parametric exci- 
tation due to single-electron effects. The terms not con- 
taining (1/2) cos(2cj m t), while having m = n, lead to 



the Duffing nonlinearity of mode 1 and 3, as described 
in the previous subsection. Single-electron mode cou- 
pling is brought about through the terms not containing 
(1/2) cos(2cj m t), while having m / n, as the amplitude of 
one mode changes the spring constant of another mode. 
The change in resonance frequency of mode m due to the 
resonance of mode n is then given by: 



1 



d 3 (N) 



dv* 



, X n )Qj n Qj ri 



(26) 



As the sign of d 3 (N)/dVg goes from positive, to neg- 
ative, to positive again, across a Coulomb peak, mode 
coupling causes softening, then stiffening, and then soft- 
ening, in the CNT spring. The change in sign from mode 
coupling across a Coulomb peak arises from the mechan- 
ically modulated average charge in the same way as the 
sign of the Duffing nonlinearity does. 



VII. ADDITIONAL EXPERIMENTAL 
OBSERVATIONS OF SINGLE-ELECTRON 
SPRING AND DAMPING EFFECTS 



In sections [IV] and [V] we showed how the resonance 
frequency and the quality factor decrease as the average 
charge on the CNT is modulated through mechanical os- 
cillation. How much the average charge changes is de- 
termined by d(N)/dV g , which in turn depends on the 
broadening IV oa d due to tunnel coupling to the leads. 
First, in subsection |VII A[ we show that, by tuning the 
tunnel rates through magnetic field, the dips in resonance 
frequency and quality factor are influenced significantly, 
and in a way that is in qualitative agreement with the 
physical picture of the damping and frequency shifts from 
our model. In subsection |VII B[ a step in the average 
charge caused by an excited state of the CNT is pre- 
sented to lead to a decrease in the resonance frequency, 
demonstrating detection of the excited state of the quan- 
tum dot using the mechanical resonator. Finally, in sub- 
section I VII C[ we demonstrate that the single electron 
spring and damping effects we observe are indeed origi- 
nating from Coulomb blockade, by studying the gate de- 
pendence of the resonance frequency in the Fabry-Perot 
conductance regime. In this regime, the charge on the 
suspended CNT segment is no longer quantized, and we 
no longer observe dips in the mechanical frequency as we 
sweep the gate. 



A. Resonance frequency shifts and damping via 
magnetic-field dependent tunnel rates 

In this subsection, we examine the quality factor and 
mechanical frequency as a function of the tunnel rates 
of the quantum dot. We tune the tunnel barriers of the 
quantum dot in a somewhat unconventional way by using 
a magnetic field parallel to the CNT. 
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The tuning of the tunnel barriers by magnetic field 
occurs through the parallel magnetic field's influence on 
the bandgap of the CNT. Figure [8^i shows the current 
flowing through the CNT as a function of gate voltage 
and magnetic field parallel to the CNT. The increase in 
the bandgap is visible as an increase in the gate voltage 
range between the hole current and the first Coulomb 
peak, which can be seen clearly by the trajectory of the 
first Coulomb peak in Fig. [HJd. Also, as the magnetic 
field is increased, the width and the magnitude of the 
Coulomb peaks decreases. This decrease in the Coulomb 
peak width is a result of an increase in the p-n junction 
tunnel barrier height by the increased bandgap at higher 
magnetic fields. 
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FIG. 8. (a) Current as a function of gate voltage and parallel 
magnetic field, showing a decrease in the height and width of 
the Coulomb peaks and a shift of the Coulomb peaks with 
increasing magnetic field for the first four charge states. The 
yellow dashed box denotes the region shown in|8]3. (b) Cur- 
rent as a function of gate voltage and parallel magnetic field, 
showing the shift of the transition from to 1 electrons with 
magnetic field (c) Cross-section of the Dirac cones of a CNT 
at a constant energy as a function of k± and k\\ , showing the 
intersection of the quantized k± with the Dirac cones, (d) 
Cross-section of a Dirac cone of a CNT, showing the change 
in the bandgap, E gap , as the parallel magnetic field alters k±. 

The increase of the bandgap with an increasing paral- 
lel magnetic field can be explained as follows^. A par- 
allel magnetic field changes the quantized wavevectors 
k± of the electrons along the circumference of the CNT 
through an Aharonov-Bohm term. The band structure 
of the CNT is determined by taking a cross-section of the 



Dirac cone at a constant k± . Figure [HJ^ shows the quan- 
tized wavevectors k± intersecting with the Dirac cones. 
Figure [8]i illustrates how the bandgap is changed through 
a parallel magnetic field. The height of the tunnel barrier 
is determined by the bandgap of the CNT, through the p- 
n junction that is formed at the interface of the metal and 
the CNT. In general, the shift of the quantization lines 
with magnetic field decreases the bandgap. This occurs 
only up until a magnetic field Bc> irac , at which point 
the quantization line crossing the Dirac point and the 
bandgap begins to increase again. From the orbital mag- 
netic moment we observe of 1.0 meV/T, together with 
the bandgap of E g = 58 meV at zero magnetic field, this 
should occur at a very large magnetic field of 29 T. Sim- 
ilar to previous reports 64 , we observe that the quantiza- 
tion line crosses the Dirac point at a magnetic field much 
smaller than expected, in this case at 0.6 T. At a mag- 
netic field above Br> irac , the p-n junction tunnel barriers 
to the quantum dot increase in height, and the tunnel 
rates to the quantum dot are significantly reduced. 

The tunnel rate of the quantum dot to the leads as 
a function of magnetic field, Ti) roa( i(B), is determined 
from the observed width of the Coulomb peaks. Figure 
shows the tunnel coupling as a function of parallel 
magnetic field, decreasing with an order of magnitude. 
To determine the tunnel coupling, a Lorentzian fit, sim- 
ilar to Eq. ([2|, was performed on the Coulomb peak at 
V b = 0.3 mV. 

In Figs. [9}d and[9]3, the resonance frequency and quality 
factor as a function of parallel magnetic field are shown, 
respectively. The resonance frequency and quality factor 
are determined by taking a frequency response at the top 
of the Coulomb peak of the charge transition from 7 to 
8 electrons, at a low bias voltage of Vb = 0.3 mV. Be- 
tween T and 9 T, the reduction of the tunnel coupling 
with a parallel magnetic field causes the resonance fre- 
quency to decrease by ~1 MHz. At a magnetic field of 9 
T, the quality factor has decreased by a factor of ~200, 
compared to zero magnetic field. 

The decrease in resonance frequency and quality factor 
with magnetic field can be explained as follows. As the 
tunnel coupling is decreased through a parallel magnetic 
field, the charge transition takes place in a smaller range 
of gate voltage, yielding a larger slope d(N)/dV g . The 
larger change in the displacement-dependent force act- 
ing on the CNT during a mechanical oscillation causes 
a larger change in spring constant and therefore in res- 
onance frequency. Also, the asymmetry of the CNT 
level between being above and below the Fermi level of 
the leads is sharper, leading to a larger retarded single- 
electron force acting on the CNT, and consequently a 
lower quality factor. 

The decrease in quality factor with magnetic field has 
to be put in contrast with the recently observed magnetic 
damping in CNT resonators^. Damping in Ref. [18] is 
measured with a CNT placed perpendicular, not parallel, 
to the magnetic field. There, Lorentz forces are the cause 
of damping, due to the perpendicular orientation of eddy 
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FIG. 9. (a) Broadening due to tunnel coupling as a function 
of magnetic field, showing a decrease in broadening with in- 
creasing magnetic field, (b) Quality factor on the Coulomb 
peak of the 7 to 8 transition at V& = 0.3 mV showing the 
quality factor decreasing by two orders of magnitude with in- 
creasing magnetic field, (c) Shift in resonance frequency on 
the Coulomb peak of the 7 to 8 transition at V& = 0.3 mV, 
showing a decrease of ~1 MHz. 

currents flowing through the CNT with respect to the 
magnetic field. In this article, damping is increased as the 
magnetic field reduces the level broadening by increasing 
the band gap, and the modulation of the average charge 
due to the mechanical motion is increased. 



B. Mechanical detection of an excited state 

In this subsection, we demonstrate that the presence of 
excited states of the quantum dot inside the bias window 
can also result in frequency shifts of the mechanical res- 
onator. Figure [TO^ l shows the differential conductance, 
dl/dVb, as a function of bias and gate voltage in the 
charge transition from to 1 electron in a stability dia- 
gram of a different device (B). An excited electronic state 
of the CNT is visible as a diagonal line inside the single- 
electron tunneling region at negative bias. Because of 
charging effects, it is not possible for both states to be 
occupied by an electron; current still takes place through 
single-electron tunneling, but now can occur through two 



channels: through both the ground state and the excited 
state. This means that the rate to tunnel onto the quan- 
tum dot, T + , has increased. When both states are al- 
lowed for tunneling, the current increases in a step, lead- 
ing to a high dl/dVb, as visible in the stability diagram. 




-4-2 2 
Gate voltage offset (mV) 

FIG. 10. (a) Stability diagram at V g = 293 mV of device 
B, showing an excited state as a diagonal line at negative 
bias voltage inside the single-electron tunneling region, (b) 
\(I — Io)\/\I — Io\max as a function of drive frequency and gate 
voltage offset at Vb = — 2 mV, showing three dips in resonance 
frequency, indicated by black arrows, with the middle dip 
corresponding to the excited state. 

The influence of the excited state on the mechanical 
motion of the CNT is examined through a measurement 
of the normalized current as a function of gate voltage 
and drive frequency at a bias voltage of -2 mV, shown 
in Fig. [To} ). Three dips in resonance frequency can be 
seen: the leftmost and rightmost correspond to the chem- 
ical potential of the ground state aligning with the Fermi 
level of the source and drain, respectively. The middle 
dip corresponds to the chemical potential of the excited 
state aligning with the Fermi level of the drain. At this 
gate voltage, the excited state causes the rate to tunnel 
onto the CNT to increase, whereas the rate to tunnel 
off the CNT remains the same. This leads to a step in 
the average charge residing on the CNT. As with the 
ground state, tunneling through the excited state causes 
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the CNT to experience a displacement-dependent force 
as it oscillates, resulting in a frequency dip. 



C. Absence of single-electron spring effects in the 
Fabry-Perot regime 

In this subsection, we examine the regime of low tun- 
nel resistance Rt for its influence on the mechanical mo- 
tion by looking at the Fabry-Perot regime 53 , for which 
Rt < e 2 /h. For quantum dots with large tunnel bar- 
riers (Rt > e 2 /h), the charge on the quantum dot 
is quantized. For sufficiently transparent barriers with 
Rt < e 2 /h, however, Coulomb blockade is destroyed 
by quantum fluctuations of the charge, and charge on 
the quantum dot is no longer quantized. In this regime, 
the conductance as a function of gate voltage still os- 
cillates due to electronic Fabry-Perot interferences, but 
the charge quantization is lost. Figure [TT|l shows Fabry- 
Perot oscillations in the measured voltage across the 
CNT as a function of gate voltage using a four termi- 
nal current bias measurement of device B. Figure [Tl^ 
displays the measured mechanical resonance frequency 
as a function of gate voltage and drive frequency. Visible 
in the plot is an increase in the resonance frequency as 
gate voltage is decreased, resulting from the electrostatic 
force from the gate that induces tension in the CNT. 
However, frequency dips do not occur at a Fabry-Perot 
oscillation. This is because the Fabry-Perot oscillations, 
in contrast to Coulomb oscillations, are not associated 
with discrete steps in the average charge. Softening due 
to the electrostatic force on the CNT quantum dot is 
therefore constant (or very slowly varying) across the en- 
tire gate range, and no frequency dips occur. 



> 

3 




voltage (V) 



FIG. 11. (top) Voltage measured across the CNT in a four 
terminal current bias configuration as a function of gate volt- 
age, showing Fabry-Perot oscillations at a bias current of -20 
nA. (bottom) Change in measured voltage as a function of 
drive frequency and gate voltage, showing the tuning of the 
mechanical resonance frequency with tension arising from the 
gate. 

resonance frequency. The occurrence of frequency dips 
in the Fabry-Perot regime is excluded due to the absence 
of steps in average charge. 



VIII. CONCLUSIONS 



To conclude, we have found quantitative agreement be- 
tween the experimental observation of single-electron ef- 
fects on the resonance frequency and quality factor of a 
CNT quantum dot resonator, and a theoretical model. 
This allows the mechanical motion of a suspended CNT 
quantum dot to be used as a probe to detect its aver- 
age charge. It is found that the mechanical resonance 
frequency and quality factor are reduced as the average 
charge changes the electrostatic force during a mechani- 
cal oscillation. At high bias, a double dip structure arises 
for both the resonance frequency and the quality fac- 
tor, which is quantitatively supported by the model. A 
model, describing single-electron induced Duffing non- 
linearity and mode coupling, leads to the finding of a 
single significant contribution. Additional experiments 
illustrating the model show that, by tuning the tunnel 
rates, the resonance frequency is reduced by ~1 MHz, 
and the quality factor by a factor of ~200. The increase 
of the average charge due to tunneling through an ex- 
cited state of the CNT also leads to a reduction of the 
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Device properties 



m 


2.6 • 10~ 21 kg 


fo 


286.82 MHz 


c 9 


2.9 aF 



Model parameters 



figure 



Abroad 
Ctot 

dC g I dx 

Qint 

N 

B\\ 
V g 

CLL 
ClR 

b L 
b R 



7 and 12 



450 GHz 
8.2 aF 
-6.5 zF/nm 

1 

T 
0.565 V 
140 GHz 
170 GHz 

288 /eV 
288 /eV 



540 GHz 
12.3 aF 
-5.2 zF/nm 
100,000 
7 
8 T 
0.9 V 
22 GHz 
22 GHz 
20 /eV 
20 /eV 



TABLE I. Properties of device A and parameters that were 
used for the quantitative fits in Figs. [4] and [6] 



Appendix A: Parameters used in the model 

For the quantitative fits describing the resonance fre- 
quency, quality factor, and nonlinear ity parameters in 
Figs. [4j [5j |6j [7| and [l2j the following parameters were 
used. 



ure 12 shows that the other contributions to the single- 
electron Duffing parameter are more than one order of 
magnitude smaller and can be neglected. We use the no- 
tation where is the contribution to a set resulting 
from the term with in Eq. 

2 
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FIG. 12. (top) Average occupation, (N), as a function of gate 
voltage, identical to the top panel of Fig. [4^,. (bottom) Other 
calculated contributions to the nonlinearity parameter, which 
are shown to be more than one order of magnitude smaller 
than a° SET . 



Appendix B: Other single-electron contributions to 
the Duffing parameter 



In subsection VIA we discussed the different single- 
electron contributions to the Duffing parameter, arriving 
at the dominance of the leading order term oP SET . Fig- 
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